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Using exact time- dependent non-singular solutions 
[Mineev-Weinstein and Dawson, Phys. Rev. E 50, R24 
(1994); Dawson and Mineev-Weinstein, Physica D 73, 373 
(1994)], we solve the Saffman- Taylor finger selection problem 
in the absence of surface tension. We show that a generic 
interface in a Hele-Shaw cell evolves to a non-linearly stable 
single uniformly advancing finger occupying one half of the 
channel width. This result contradicts the generally accepted 
belief that surface tension is indispensable for the selection of 
the i-width finger. 



The problem of the finger width selection was posed in 
1958 by Saffman and Taylor [Q in their study of displace- 
ment of oil by water in a Hele-Shaw cell. This cell consists 
of two parallel glass plates separated by a thin gap occu- 
pied by a viscous liquid which is pushed by a less viscous 
one. This simple device is useful for modeling flows in 
porous media, the study of which is vitally important for 
many applications. Flows in uniform porous media and 
in the Hele-Shaw cell are described by the same Darcy 
law: v = — Vp, where v is fluid velocity, and p is pres- 
sure. For the Hele-Shaw cell, this equation follows from 
the Stokes equation averaged over the direction perpen- 
dicular to the parallel plates. 

Saffman and Taylor Q observed that an almost planar 
initial oil/water interface becomes unstable and gives rise 
to many competing fingers which eventually evolve to a 
single uniformly advancing finger occupying one half of 
the channel width, if the surface tension is very small. 
But, as was analytically shown in the same paper |Q, 
any finger width is possible. So the selection problem was 
stated: why does Nature choose the width of one-half? 

This problem appeared to be universal, i.e. the same 
selection phenomenon is common for displacement of var- 
ious viscous liquids by less viscous ones for immiscible 
incompressible liquids. This problem is related to the 
problem of pattern selection in noncquilibrium phenom- 
ena, which has been of much subsequent interest . 

When the viscosity of water is negligible compared 
with oil viscosity, the mathematical formulation of this 
problem in the absence of surface tension has the form: 



y-p — (in the oil domain) , 

p = —x if x — > +oo (oil pushed to the right), 

p = (at the oil/water interface), 

d n p = for y = ±7r (at the channel walls), 

V n = —d n p (at the oil/water interface), 



(1) 



where V n is the normal velocity of the interface, d n p is 
the normal component of the pressure gradient, and the 
channel width is chosen to be 2n in our scaled units. 

The solution of the system (Q) , describing a finger mov- 
ing in the x-direction with velocity 1 /A and occupying the 
portion A of the Hele-Shaw channel width, is jy 



= 2(1 - A) log[cos(y/2A)] + i/A . 



(2) 



We parametrize the moving interface z(t, <j>) = x{t, <£>) + 
iy{t, <j>) at time t by the parameter cf>e [0, 2tt]. After the 
shift y — > y — 7rA, Eq.(2) can be rewritten as 



z(t, <f>)=t/\+(l- a)i<j) + a log[(e^ - 1) /2< 



(3) 



where A = 1 — a/2. The system (1) can be reduced to 
what we call the Laplacian growth equation (LGE) for 
the moving front z(t, <ft) (see || and references therein): 



Jm(ztz^) 



(4) 



Here the bar denotes complex conjugate, zt and z^ are 
partial derivatives, and the map z(t,<j>) is conformal for 
Im cj) < 0. One can easily see that z(t, <f) given by (3) 
is the traveling- wave solution of the LGE given by (0). 
The finger width A is here a free parameter, while exper- 
imentally it is always 4. What determines A? 

In jj| Saffman and Taylor proposed that surface ten- 
sion between the two fluids would solve the selection 
problem. Since then, it has been widely accepted that 
the inclusion of surface tension is the only way to se- 
lect the most stable finger width, and much work was 
done toward solving the selection problem in this way 
(see books H and references therein). While mathemat- 
ically non-trivial and challenging [Q, this activity, espe- 
cially intensive in the 1980s, was nonetheless successfully 
completed 0j and summarized in H. In short, several 
groups in 1986-1987 0] confirmed (using expansion "be- 
yond all orders" and reduction to a nonlinear eigenvalue 
problem) numerical evidence || of the discrete spectrum 
of A, decreasing to 1/2, in the limit of low surface ten- 
sion. Surface tension was claimed to be responsible for 
the selection: equating it to zero would make this analy- 
sis senseless. 

We see two reasons for the absence of attempts to ex- 
plain the selection without the inclusion of surface ten- 
sion. First, because of the absence of analytic time- 
dependent solutions, all selection studies were focused 
on linear stability analysis of a steady-state traveling fin- 
ger in the presence of surface tension, which is the main 
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physical factor neglected when the continuous family of 
fingers (2) was derived JjJ. The second reason to include 
surface tension stems from the observation that almost 
all exact zero surface tension solutions of this problem 
obtained before 1994 B exhibit finite-time singularities 
(cusps). Due to the belief that these solutions are gen- 
eral and thus capture main features of this problem, it 
was concluded that to reach long times, it is necessary to 
include surface tension to eliminate singularities M . 

In 1994, we reported jj| a new class of exact time- 
dependent solutions of the LGE (|4|) having the form: 

N 

z{t, <f) = r(t) + ificj) + J> fc lo S( e# ~ **(*)) . (5) 



fc=i 



where jj, = 1 — Y2k=i ak ' ak = const, and |a&| < 1. With 
some (quite modest) constraints on {ak}, these sohitions 
remain non-singular and analytic for all times (no cusps) 
pL The time dependence of a.k(t) and r(t) is given by 



JV 



0k = z{t,i loga fc ) = r - (1 - ^a ; )loga fc (6) 

i=i 

N 1 
+ > cei log(- ai) = const, and 

1=1 

1 N 1 N 

t + C= (l--^a fc )r+-^a fc log(a fc ), (7) 



k=l 



k=l 



where fc = 1,2, ...,JV and C is a constant in time [|10| . 
Eqs. (6) and (7) follow from the substitution of (5) into 
(4). 

All ak are located inside the unit circle and, if the same 
holds for the roots of z$, then z(t,(j>) is conformal for 
Im (j> < 0. We called these solutions TV-finger solutions, 
since they describe the evolution of N fingers. The class 
of solutions (5) contains all previously known exact solu- 
tions H , [jll| , |Q , including those which diverge in finite 
time. The subclass of these solutions without finite-time 
singularities, was also shown to be dense in the space of 
all analytic curves |l5| ] . The dynamics of an arbitrary ini- 
tial interface can be faithfully described within this class. 
In addition to possessing these attractive mathematical 
properties, these solutions are very physical: they de- 
scribe tip-splitting, side-branching, competition, coarsen- 
ing and screening of growing fingers which are observed 
in all known experiments and simulations. 

The following geometrical interpretation of the con- 
stants {afc} and {{3k} @ is of great help: the complex 
number (3k — otk log 2 is the location of the fc-th stag- 
nation point which the interface does not cross (tips of 
white grooves in Fig. [l] and Fig. [2]), but approaches expo- 
nentially slowly, namely proportional to exp(— r/Re ak)- 
Near the fc-th stagnation point, a groove with parallel 
walls originates, with width 7r|a:fc| and angle with respect 



to the horizontal axis Ok = arg ak- In terms of these 
stagnation points and grooves given by constants {Pk} 
and {ak}, all the dynamical features mentioned above 
are especially clear (see Fig. |l|). These grooves merge 
and finally coalesce to a single growing finger in accor- 
dance with all known experiments and simulations KLpJ . 
Formally this means that a generic initial interface given 
by the iV-finger solution (5) necessarily evolves to a single 
uniformly advancing finger (see [Q for details). 
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FIG. 1. Geometrical interpretation of the complex con- 
stants of motion <%k and k = 1, N. 

In this paper, we will solve the selection problem ana- 
lytically, starting from first-principles, in the absence of 
surface tension (while interfacial tension is required to 
linearly stabilize the finger, as previously shown ) ■ 

It is known |ll| that the development of a single finger 
with a relative width A = 1 — a/2 in the long-time limit 
is described by 



z(t, 4>) = r(t) +fj,i(f> + a log (e 10 - a(t)) 



(8) 



> OO , T — 

a)] , where c 



where 0<a<2, < a < 1. For t 
2b/ {2 - a) and a(t) = 1 - cexp[-2i/a(2 
and a are constants in time. Choosing the width of the 
Hele-Shaw cell to be 27r, we have z[t, 27r) — z(t, 0) = 2iri, 
because of the periodic boundary conditions. Calculating 
z(t, 2ir) — z(t, 0) from (||) and using the fact that \a\ < 1, 
we obtain that 2iri = 2 j ni[i-\-2'Kia = 2ni(iJ, + a), or finally 
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(9) 



Then we note that the second term in the right-hand side 
of (g), namely \x%4>, is the limiting value of the logarithm 
with a logarithmic pole e located at zero: 



u,i<j) — /i lim log (e 



itp 



(10) 



Let us perturb the interface (8), corresponding to £ = 
0, by the initially small non-zero e, < e <C 1. The 
perturbed interface will have the form 

£(t,<£)=r + Mlog(e^-e)+alog(e^-a) , (11) 

Now one can easily see that the value e = (and thus 
the finger described by (8) with \i ^ 0) is unstable. The 



2 



point is that Eq. ([□]) is exactly the N = 2 case of the N- 
finger solution (5) of the LGE (4). As one can see from 
(6) and (7), a(t) and e(t) merge at unity when t — * oo: 

a=l-he- T , e = l-l 2 e- T , (12) 

where the constants l\ and I2 are determined by 



JV 



fa = fj,log(h + l 2 ) + a\og(2h) , 
f3 2 = /ilog(2Z 2 ) +a\og{h +l 2 ) ■ 



(13) 



In view of (|12|), we substitute 1 for both a and e with the 
accuracy 0(e _T ), and thus obtain from (11), for t — > 00 



z(i,0)=T+Ou + a)log(e^-l) 



(14) 



Let us interpret the result (|l4j) : due to the instability of 
the initial finger with width A = 1 — a/2, the new 
finger described by (Oh has been formed. Its width is 



X" 



l-( M + a)/2 = l/2 i 



(15) 



in accordance with the condition (|9|) that fx + a = 1 . 

Let us perturb the finger (||) in a more general way, 
than we did in (11). We note that 



Liitf) = lim ^ 4 log (e^ - e k ) 

all e t ->0 *-~> V 7 



if y\ —1 Sk = fi • Choosing all to be nonzero, we 
rewrite the finger (||) in a perturbed way as 

JV 

z(t, <f>) = t + S k log (e* - e k ) + a log (e* - a) . (16) 
fc=i 

Equation ( [Tif ) is the (jV+l)-finger solution (5) of the LGE 
(4) with dynamical conditions (6) and (7). Because of the 
density of the subclass of smooth solutions given by (5) 
we conclude that (|l^) describes a general perturbation of 
the finger (||), if N is large enough. As follows from (6) 
and (7), generally all logarithmic poles in the absence of 
finite-time singularities merge [similarly to (12) for N — 
2] in the long time limit to 1 with exponential accuracy 
0{e- T ) @, where 



, JV 



fc=l 



(17) 



Because of this merging near the unit circle we substitute 
1 for all tk{t) and a{t) in (Hq) when t — > 00 and obtain 



(18) 



JV 

z{t,</>)=T+(j28 k + a) log(e^-l) 



fe=i 



This formula describes the single finger formed from (8) 
under the perturbation (|l6|). Its width is 



= 1 - E 4 + 7 2 



(19) 



k=l 



which is exactly one half since 5^ = /j, (see above) 

and fj, + a = 1 by an argument analogous to (m , so that 



A now = 1 - (ji + a)/2 = 1/2 



(20) 



Both ( |20D and (g_5|) indicate that, for obtaining instability 
of non-i-width finger and formation of the ^-width finger 
in a long-time limit, surface tension is not needed. 

Now we start from an arbitrary initial interface in 
terms of non-singular solutions (5), (Re otk > 0). Be- 
cause of the coalescence of all initially non-zero poles at 
at the unit circle, the long-time limit of (5) is given by 

a finger (3) with a width of (2 - J^k=i "fe)/ 2 - The N ~ 
finger solution (5) is the limit of (N + A:)-finger solution 
expressed by the same equation, but without the term 
This limit corresponds to equating k of the poles to 
zero, and this value of zero can be easily shown (in the 
same way as above) to be unstable for all of these k poles. 
Because of the density of these solutions in the class of 
all analytic curves (l|, this (N + /c)-finger solution can 
be arbitrarily close to any analytic interface. Again, all 
a k (t) merge to 1 in the limit t —> 00 [as stated earlier and 
proven in J14|]. Thus we have in the long-time limit 



JV+Jc 



l) = T + log I 



a i<t> 



1 



k=l 



Here we used ~}2ik=i ak = 1> smce 27ri = z(t, 2ir)~z(t, 0). 
So, we have demonstrated that initial interfaces evolving 
to the non-i-width finger are unstable with respect to 
formation of the ^-width finger, which thus is shown to 
be the only attractor for all generic moving fronts in the 
Hele-Shaw cell represented by (5). (Solutions with sev- 
eral parallel fingers forming in asymptotics can also be 
easily shown to be unstable). The dynamics of the tran- 
sition from an arbitrary interface to the ^-width finger is 
exactly described by the set of transcendental equations 
(6) and (7) which involve only elementary functions. This 
selection of one half is in agreement with known experi- 
ments and simulations in the limit of low surface tension. 

Now we will extend these results obtained for peri- 
odic boundary conditions to the more physical "no-flux" 
boundary conditions (no flow across the lateral bound- 
aries of the channel) . This requires that the moving in- 
terface orthogonally intersects the walls of the Hele-Shaw 
cell. However, unlike the case of periodic boundary con- 
ditions, the end points at the two boundaries do not nec- 
essarily have the same horizontal coordinate. This is also 
a periodic problem where the period equals twice the 
width of the Hele-Shaw cell. The analysis is the same as 
before, but now only half of the strip should be consid- 
ered as the physical Hele-Shaw channel, while the second 
half is the unphysical mirror image (see Fig. |^) . 
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To be brief we will perform the analysis for N = 2, 
but one can trivially extend it to an arbitrary N. The 
easily obtainable extension of (@) for the development of 
a single finger with the width A = 1 — (a.\ + a 2 )/2 is 

z(t, 4>) = r(t) + flic/) + ax log (e 10 - Oi(i)) 
+a 2 log(e i * + a a (t)) , (21) 

where < ai,2 (t) < 1. The generalization of (9) is 

H = l-{a l + a 2 ). (22) 



A.y 




FIG. 2. Only half of the periodic cell should be considered 
as the physical Hele-Shaw channel, while the second half is the 
unphysical mirror image. 

We note that 

Hi<j>=(l-ai + 6) lim log {e*+ - e x ) (23) 

i ei— +0+ 

+(- -a 2 -8) lim log (e^ + e 2 ) . 

Substituting this into and allowing ei j2 to be func- 
tions of time, we note that z(t, <j>) is the N = 4 solution 
(5) of Eq.(4), where 8 and ei )2 must be real. If |£| < 1/2, 
all four logarithmic singularities ai j2 (i) and ei, 2 (i) merge 
to 1. So we have in the long-time limit 

z(t,0)=T(i) + (i + £)log(e^-l) (24) 

+(±_*)log(e* + l) . 

This describes a finger moving between two grooves with 
widths 7r(l + 26) /2 and 7r(l — 28)/ 2 respectively (see the 
geometrical interpretation above). Thus the portion of 
the channel width occupied by the moving finger is 

= 2tt - tt(1 + 2^/2 - tt(1 - 2^/2 = 1 

2tt 2 ' 1 j 

as before. So, for the no-flux boundary conditions, we 
have obtained the ^-width finger as expected. In exper- 
iments, the finger in the long-time limit is centralized in 



sense that axes of symmetry of the finger and Hele-Shaw 
cell coincide. This corresponds to the condition (5 = 0. 

The finger with a width of one half is nonlinearly 
stable with respect to generic perturbation of logarith- 
mic type (5). Namely, the shape of the finger can be 
destroyed at the initial (linear) stage, but eventually the 
i-width finger will be restored, because of the coalescence 
described above. (Of course, the i-width finger is lin- 
early unstable without interfacial tension in accordance 
with previous studies 0||,|7]].) 

Regarding surface tension, we think that while math- 
ematically still singular (because a small number multi- 
plies the highest derivative) , physically surface tension is 
a regular perturbation for this problem (unless very high 
curvatures exist which surface tension suppresses). 

In conclusion, we have analytically solved the finger 
selection problem in the absence of surface tension. By 
using the non-singular exact solutions of the LGE (4), 
we have demonstrated that the i-width finger is the only 
attractor for all generic moving fronts in a Hele-Shaw cell. 

The author thanks H. Makaruk, W. MacEvoy, J. Pear- 
son, V. Lvov, I. Procaccia, M. Feigenbaum, D. H. Sharp 
for helpful discussions. 
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Note added after publication: In the recent paper by 
A. P. Aldushin and B. J. Matkowsky (to be published 
in Appl. Math. Lett.) the selection of A = 1/2 in the 
absence of surface tension was shown to be consistent 
with minimum of assumed functional for the steady-state 
solution. 
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